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Abstract We present a new automated algorithm to identify, track, and charac- 
terize small-scale brightening associated with solar eruptive phenomena observed 
in Ha. The temporal spatially-localized changes in chromospheric intensities can 
be separated into two categories: flare ribbons and sequential chromospheric 
brightenings (SCBs). Within each category of brightening we determine the 
smallest resolvable locus of pixels, a kernel, and track the temporal evolution of 
the position and intensity of each kernel. This tracking is accomplished by iso- 
lating the eruptive features, identifying kernels, and linking detections between 
frames into trajectories of kernels. We fully characterize the evolving intensity 
and morphology of the flare ribbons by observing the tracked flare kernels in 
aggregate. With the location of SCB and flare kernels identified, they can easily 
be overlaid on top of complementary data sets to extract Doppler velocities and 
magnetic field intensities underlying the kernels. This algorithm is adaptable to 
any dataset to identify and track solar features. 

Keywords: Flares, Dynamics; Chromosphere, Active 
1. Introduction 

Active regions on the Sun can produce a variety of distinct dynamic features 
including sunspots, flares, prominences and occasionally mass ejections; as such 
solar active regions are a major source of adverse space weather. Historically, 
these active regions have been identified and characterized manually. To better 
understand the dynamics of solar active regions and their associated phenom- 
ena, data from different sources must be utilized and a large dataset consisting 
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of numerous active regions must be analyzed. This necessitates an automated 
approach to correlate several seemingly disparate data sets. 

In the past, a few different automated methods have been developed to si- 



multaneously analyze flare ribbons across multiple sources of data ( 


Qu et al., 


2003 


Maurya and Ambastha, 2010 


Gill, Fletcher, and Marshall, 2010 


), however 



sequential chromospheric brightenings (SCB) have not been studied in such an 
automated fashion. First observed in 2005, sequential chromospheric brighten- 
ings (SCBs) appear as a series of spatially separated points that brighten in 



sequence ( Balasubramaniam et al., 2005). The sequential nature of the point 



brightening gives the appearance of a progressive traveling disturbance. SCBs 
are observed as both single and multiple trains of brightenings in association with 
a large-scale eruption in the chromosphere or corona. The loci of brightenings 
are seen in a time-series of full-disk Ha images as emerging predominantly along 
the axis of the flare ribbons. Physically, SCBs are closely correlated with solar 
flares, coronal restructuring of magnetic fields, halo CMEs, EIT waves, and 



chromospheric sympathetic flaring ( jBalasubramaniam et al., 2005 ). Pevtsov, 



Balasubramaniam, and Hock (2007) show that these physical correlations are 
consistent with several properties of chromospheric evaporation. 

This paper presents a new automated method of identifying and tracking 
SCBs and associated flare ribbons. This tracking technique is different than 
previous flare tracking algorithms in that it identifies and tracks subsections of 
the flare from pre-flare through the impulsive brightening and the exponential 
decline. This allows measurement of the temporal variation in intensity and 
position, which can then be combined with other data to infer the changes 
occurring to the underlying photospheric magnetic field as well as derived Ha 
Doppler velocity of each individual flare kernel. This tracking algorithm is also 
adapted to identify and track the temporal evolution of the ephemeral SCBs 
associated with the flaring region. Section [2] describes the data sets used to train 
the algorithm and the image processing involved before the detection routine can 
be applied. Section [3] explains the method of detection used to identify and track 
the flare and SCB kernels. Section [4] presents the capabilities of applying this 
detection algorithm to flaring events. Finally, Section [5] discusses the implications 
of these results and the future direction of this project. 



2. Data 

This project utilizes Ha (6563 A) images from the Improved Solar Observing 
Optical Network (ISOON) telescope to train the algorithm in identifying and 
tracking flare ribbons and SCBs. ISOON produces a 2048 by 2048 pixel full- 



disk image with a one-minute cadence (Neidig et al, 19981. These full-disk 
chromospheric images have a 1.1 arc-second resolution and are normalized for 
quiet Sun intensity and corrected for atmospheric refraction (Figures [T]A., [5]A., 
and|3]A.). At the same resolution and cadence, ±0.4 A off-band Doppler images 
are also recorded. 

In this study, we chos e three flares (Figures [l| [2{ and [3]) with previously 



identified SCB events by Balasubramaniam et al. (2006) as a sample set to 
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Table 1. The events selected by this study to analyze. Each 
event was previously identified as having SCBs associated with 
the flare and had a two-ribbon topology. 



Date 


Start Time (UT) 


Flare Class 


CME 


2004 November 9 


16:30 


M 


Halo 


2005 May 5 


16:03 


C 


Halo 


2005 May 13 


16:33 


M 


Halo 



train the detection algorithm (see Table [I]). Each of the events selected had 
halo CMEs associated with the flare, a filament eruption, and had a rough 
two-ribbon topology. Both 6 and 13 May 2005 events were located near disk 
center while the 9 November 2004 event was near the eastern limb. To achieve 
complete coverage of the flare's eruption, images from ±3.5 hours from each 
eruption's start time are extracted from the archive. This yields an image cube 
of ^400 images for an individual event. To characterize the flare eruption, three 
methods were used: a relative intensity curve, a maximum intensity curve, and 
the Geostationary Operational Environmental Satellite (GOES) hard and soft 
x-ray fluxes (Figures [l]C, [2]C, and|3]C). The relative intensity curve is a ratio of 
the maximum to mean intensity for each of the Ha images: 

j ^max ^mcati \ 

= 1 +7 ' 1 ' 

1 max ~ -'mean 

where I r is the relative intensity, I max is the maximum intensity, and J meall is the 



mean intensity of each Ha image (McAteer et al., 2002 McAteer et al., 2003). 



The maximum intensity curve is / ma x for each image. The GOES Soft (1.0 — 8.0 
A) and Hard (0.5 - 4.0 A) curves are x-ray intensity curves from the GOES 14 
satellite. 

Before the detection analysis is applied, each image is preprocessed (see Fig- 
ure]^]). Preprocessing of line-center Ha images begins by removing limb darkening 
from the calibrated full disk ISOON image. This is accomplished by subtracting 
a standard limb darkening profile determined by the spherical geometry of the 



solar surface (Hestroffer and Magnan, 1998 Neckel and Labs, 1994). Each image 



is then de-projected into conformal cylindrical coordinates which removes the 
projection effects of imaging the solar sphere. As a consequence, latitudes above 
85 degrees north or south are significantly distorted and are cropped off. In this 
cylindrical projection, the entire solar limb is distorted. Any region of interest at 
the limb is significantly distorted as well leading to added noise in the analysis; 
for an example of an event on the limb, see Figure [3] The de-projected image 
is then normalized to the mean background level of the quiet Sun so that any 
flare brightening is comparable between images. Next, the full-disk de-projected 
image is cropped down to isolate the flaring region (Figures [l}3, [2^3, and|3^). 
The flaring regions in each image are aligned using a cross-correlation algorithm. 
This eliminates the rotation effects of the Sun so the flare can be observed as 
if it was stationary on the solar surface. An algorithm to de-stretch each image 
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is applied to the image cube to co-align subsections of the image. The mean 
displacement of any given subsection due to de-stretching over time is nearly 
zero. A net zero displacement implies the algorithm is not removing any proper 
motion observed in the time series. A discussion of the methodology behind the 
de-stretching algorithm can be found in November and Simon ( 1988 ). The entire 
image preprocessing of a flare is automated up to this point, given a library of 
ISOON images and the coordinates of the active region that is of interest. The 
last step in preprocessing is removing bad frames. Images that contain peculiar 
observations due to clouds or other factors are discarded interactively. 



3. Tracking Algorithm 



When viewing a sequential time series of images covering an erupting flare, 
several physical characteristics of the evolving ribbons are immediately apparent: 
the ribbons separate, brighten, and change their morphology. Along with the 
erupting flare, SCBs can be observed brightening and dimming in the vicinity of 
the ribbons. This section describes the techniques and methods used to extract 
quantities of interest such as location, velocity, and intensity of the flare ribbons 
and SCBs. Section [3TT] describes the technique used to identify and track bright 
kernels in a set of images. In this context, we are defining a kernel to be a small 
locus of pixels that are associated with each other through increased intensity 
as compared with the immediately surrounding pixels. Each kernel has a local 
maximum, must be separated from another kernel by at least one pixel, and does 
not have any predetermined size or shape. Section [3.2| describes some of the image 
processing required to extract characteristics of flare kernels and Section |3.3| 
describes the technique for SCB kernels. 



3.1. General Kernel Detection and Tracking 



The algorithm developed to identify, track and extract physical quantities in the 
ISOON dataset does so in seven logical steps: thresholding, feature isolation, fea- 
ture enhancement, locating the position of candidate kernels, eliminating "false" 
kernels, linking time-resolved kernels together into trajectories, and extracting 
physical quantities by overlaying the trajectories over other datasets (numbers 
3 - 9 in Figure To aid in this detection, tracking software developed by 
Crocker, Crier and Weeks (Crocker and Grier, 1996) was used as a foundation 
and modified to fit the needs of this project. Crocker's software is available online 
at www. physics.emory.edu/~weeks/idl/ 

Analysis begins with thresholding the preprocessed Ha images. A threshold 
is taken of the background-normalized images eliminating any features that 
are dimmer than a specified intensity. Because the mean intensity of the flare 
ribbons are significantly brighter than SCBs, the specific thresholding level is 
different for each feature extraction. Section [372] describes the process for flares 
and Section l3~3l for SCBs. 

Next, to isolate features a spatial bandpass filter is applied to suppress pixel 
noise with a characteristic length scale of a pixel. Noise reduction is necessary to 
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remove the smallest changes in brightening while preserving the general structure 
of the images. Selecting the length scale for the bandpass filter requires prior 
knowledge of the features' characteristics. With ISOON images, flare ribbons 
appear as ~10 pixels wide and ^100 pixels long - significantly larger than the 
bandpass filter. SCBs appear characteristically smaller in ISOON images with 
diameters of a few (3 — 8) pixels. The noise filter eliminates some of the smallest 
SCBs but is a worthwhile tradeoff for reduced noise levels. 

The process of feature enhancement is different whether looking at SCBs or 
flare ribbons. In either case, it is a process of folding a processed image back into 
the original to enhance either the flare ribbon or the SCB. Section [3~2| describes 
the process for flare tracking and Section [373] for SCBs. 



3.1.1. General Kernel Position Determination 

Determining the location of kernels is an iterative process. The corresponding 
flowchart is shown in Figure [5] The process begins by identifying local intensity 
maxima within the enhanced image, I', as candidate kernels. A pixel is consid- 
ered a candidate centroid of a kernel if a brighter pixel does not exist within a 
radius of ui pixels of the candidate. In general, lj is the characteristic radius of the 
kernel the algorithm is looking for. We required an lj to maximize the number 
of detections along a flare ribbon as well as the number of frames in which 
a detection could be tracked. Without specific knowledge of what size a flare 
kernel or SCB kernel could be, lj was iterated until both of the conditions were 
met. To maximize the number of detection regions and minimize the number of 
untraceable kernels, both quantities were plotted versus mask radius for three 
flares used in the training set. A mask with an 8 pixel radius (lj = 8) was 
consistently found to be the minimum radius where the diagnostic curves crossed. 

Next the brightness- weighted centroid of pixels, (x, y), for each candidate 
kernel is calculated with the local maxima having pixel coordinates (x, y) and 
intensity A(x, y): 

l)=— E (l)A{x + i,y + j). (2) 

mo is the integrated intensity of the candidate kernel: 

m = ^2 A(x + i,y + j) (3) 

t 2 +j 2 <u 2 



and oj is the radius of the mask which was defined above (Crocker and Crier, 



1996). Improper background subtraction can bias the centroid toward the center 
of the mask and away from the kernel's brightness center by adding weight to 
pixels that are not part of the features being identified. 

An offset between the brightness weighted centroid and the local maxima 
is then defined by subtracting the position of the brightness centroid from the 
position of the local intensity maximum: 

;)-(;)-(;)■ 
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Next, a refined centroid of the candidate kernel is calculated by: 

(x', y ') = (s+|,y+|). (5) 

This refined centroid is located halfway between the local intensity maximum 
and the brightness-weighted centroid. 

To evaluate the quality of the candidate kernel, the discrepancy between the 
local maximum and centroid is examined. If either |e x | or |ey| (the offset as 
described by Equation [4]) is greater than 0.5 pixels from the local maximum, 
the filtering mask is shifted by 0.5 pixels in the direction of (e x , e y ) and the new 
kernel is characterized. This iteration has the effect of guaranteeing that the 
local maximum and the centroid of the kernel are within half a pixel of each 
other. 

The final step in determining the kernel's position is to eliminate unwanted 
kernel detections or "false kernels". A threshold is applied to the integrated 
intensity of each candidate kernel. This intensity filter is intended to eliminate 
large, dim detections that may have a local maximum but are not part of the 
features in which we are interested. The threshold for this filter was determined 
empirically to eliminate marginal detections while retaining the bulk of the 
detections. Each detection now has position (x^y'j) and an integrated intensity 
rno.i t° characterize the kernel. The results of this detection process can be seen 
in the lower panels of Figure [6] 

3.1.2. General Trajectory Linking 

After the positions of all kernels within the series of images have been located, 
the next step is to associate kernels between frames to form trajectories. This 
presents a problem in determining which kernel in image I'„+i is most likely 
to correspond to a given kernel in the preceding image, I' n . It is complicated 
by the number of particles that are being tracked simultaneously. We simplify 
the problem by requiring that each kernel be identified with at most one kernel 
in the previous image. Since kernels are nominally indistinguishable from each 
other, a kernel pair's physical proximity between frames is a reasonable way to 
establish a correlation between the two. The following algorithm for trajectory 
linking is motivated by the dynamics of noninteracting Brownian particles; see 



Crocker and Grier ( 1996 1 for a more complete discussion. 

The probability that a single particle with Brownian motion will diffuse a 
distance 6 in time r is 

P(8W) = — 1=— exp (~^-) (6) 



where T> is the self-diffusion coefficient of each particle ( Crocker and Grier, 1996 ). 
For a system of N noninteracting particles, the diffusion probability distribution 
is the product of the distribution for each single particle: 

rm\T)=(—) -pl-V--). (7) 




\A-kVt 



SOLA: SCB_Techniques.tex; 8 August 2011; 1:29; p. 6 



Tracking SCBs and Associated Flare Ribbons 



Thus the best a ssignment of trajectories is one that maximizes P({Si}\r) or 
minimizes &i (Crocker and Grier, 1996). While merging or converging kernels 



result in a loss of tracking memory, we are able to visually discern any such 
occurrences by animating the trajectories. 

Calculating P({Si}\r) for every possible combination of trajectories would 
need 0(N\) computations which is unfeasible and not necessary in this case. We 
know that a flare or SCB is going to stay fairly localized and not jump sporadi- 
cally from one location to another. Thus we can define a characteristic diffusion 
length L to be the maximum distance a kernel can travel between frames. This 



diffusion length is different between flares and SCBs (see Sections 3.2 and 3.3 
for the specific differences). With a small enough diffusion length, most kernels 
have only one possible association between two frames. 

If all kernels remained present throughout the series of images, then associ- 
ating trajectories would be significantly simpler. Unfortunately this is not the 
case at all since we are interested in tracking the flare from its inception through 
decay. To account for a kernel disappearing and reappearing, "missing" links 
are assigned to a kernel that does not have a counterpart in the next frame. A 
diffusion distance of <5 = L is assigned for the purposes of calculating P({Si}\r). 
The last known location of a "missing" particle is retained for a predetermined 
number of frames. That number of frames is different for flares and SCBs (see 



Sections 3.2 and 3.3 for the specific differences). 

Linking detections into trajectories is only possible if the diffusion distance 
6 is significantly smaller than the typical spacing between kernels a. If 6 is 
greater than ~ | then trajectories are inextricably confused between images. 
This property of confusion provides a upper bound for kernel motion which 
provides insights into physical properties (see Section [4]). 

Linking the initial kernel detections into trajectories provides another op- 
portunity to filter out inconsistent or unsustained detections. A kernel lasting 
only one or two frames is a weak detection and can be filtered out. A temporal 
trajectory filter eliminates fleeting bright points that appear in a few frames but 
are not associated with feature being tracked. The specifics of this filter appears 
in Sections 13.21 and EP1 

The last step in the kernel identification and tracking algorithm is to extract 
physical properties of the data. To accomplish this extraction, the locations and 
parameters of the kernel trajectories are placed over the preprocessed dataset. 
Since the original Ha images were modified to aid in the detection process in 
thresholding and feature enhancement, the tracked kernels are placed over the 
un-enhanced Ha images to extract the intensity underneath the each kernel. 
Readers who arc interested in applying these methods are welcome to contact 
the lead author. 

3.2. Flare Ribbon Extraction 

Within the general feature detection and tracking framework presented above 
and outlined in Figure [4j there are several situations where the algorithm is 
customized to the features that it is trying to extract. This subsection will explain 
the specifics for extracting flare ribbons from the Ha dataset. 
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The first flare-specific customization is applying a threshold level of 1.35 to the 
preprocessed images where quiet Sun disk center intensity is set to unity. This 
eliminates any features that are approximately less than one standard deviation 
above background intensity. This corresponds to keeping features with a 35% 
increase in brightening over the background intensity and since all images are 
normalized to their background intensity the threshold value holds true for any 
ISOON data set. 

Next, the flares are isolated and enhanced. To isolate flares, a Laplacian 
operator is applied to the thresholded images. In this application, the Laplacian 
operator acts as an edge enhancement algorithm making smooth transitions on 
the image disappear while rapid changes in intensity stand out. The images 
operated upon are then added back into the original image, 

I' =1+ V 2 Z, (8) 

where I is the image after thresholding and the bandpass filter, and X' is the 
resultant image used for the detection of flare kernels (Figure |6| . By adding the 
Laplacian image back to the original, the edges and peaks in the flare ribbons 
are enhanced while the physical dimensions of the active regions stay the same. 

When linking flare kernels into trajectories, we define a characteristic diffusion 
length L as the maximum distance a kernel can travel between frames. For 
flare kernels, L = 4 pixels. This diffusion length was empirically found to be 
the maximum distance a flare kernel could travel to minimize the confusion of 
tags when associating kernels between images. The typical spacing between flare 
kernels, a, is approximately 8 pixels, thus L w |. 

Lastly a filter is applied to remove unwanted trajectories by evaluating their 
temporal integrity. The previous 40 frames are considered when linking flare 
kernels into trajectories. This means that a kernel can go "missing" for 40 images 
and be associated together. Using the fact that flare ribbons are relatively long- 
lived phenomena, lasting more than an hour, a filter is applied to eliminate short- 
lived detections. Any trajectories that last less than 20 frames are eliminated. 
This filter eliminates several off-ribbon detections which are associated with the 
eruption but do not characterize the evolution of the flare ribbons. 

3.3. SCB Extraction 

SCB detections are completely eliminated in the flare detection and tracking 
process. A separate process is required to extract SCB kernels and characterize 
this different type of brightening that occurs during a flare eruption. This sub- 
section details the customizations to the detection algorithm needed to extract 
SCBs. 

To begin to isolate the SCBs from the original background normalized images, 
a threshold of 1.2 is taken. This threshold eliminates all features except where 
intensities greater than 20% above background level. 

To isolate the SCB features in a series of Ha images, a running difference 
image series is created. In general a running difference image series eliminates 
long-lived features while enhancing rapidly changing intensities. In this case, 
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the current image in the series is a binary mask of all the pixels above the 
1.20 threshold. The binary image is necessary to force all brightening above the 
threshold to a uniform value and includes both SCBs as well as flare ribbons. 
Next, a median image is defined by taking a median of the previous 11 frames. 
The running difference image is calculated by subtracting the median image from 
the current binary image, 

71 -Z*rt -^ti— ll:n— 1; (9) 

where I is the binary image, I is the median image, and I' is the running 
difference image. The 11 minute duration of the median image was experimen- 
tally determined to best maximize the number of SCBs detected over the three 
training flares examined. The results of this process expose the locations of the 
SCBs while subtracting off the longer-lived flare ribbons (Figure [7]) . 

The isolated SCBs are then enhanced through a morphological transform. The 
transform performs an erosion operation followed by a dilation operation with 
a 3x3 pixel plus-shaped structuring element. This transform has the effect of 
removing small regions of noise in the running difference image while preserving 
the size and shape of the features in the image. 

Using the binary running difference images, the SCB's local maxima are 
identified as outlined in Section [3.1.1| However since the images being analyzed 
are binary images, the local maximum is coincident with the centroid of the 
bright points within the mask which significantly decreases processing time. The 
kernels are then run through the trajectory linking algorithm with a diffusion 
length L = 10. This diffusion length allows the centroid of the SCB kernel to 
move no more than 10 pixels between frames. 

When correlating SCB kernels between frames into trajectories, only the 
previous image is searched. A temporal filter is also applied to the trajectories 
requiring at least two frames. What remains are kernels that last at least two 
minutes and do not "go missing" at anytime. Lastly, a spatial filter is applied 
to eliminate any SCB kernels that are identified as within the "flaring region." 
This region is defined by co-adding 80 of the thresholded images used in the 
identification of flare kernels centered around the peak of the flare intensity. All 
non-zero pixels are nagged as the "flaring region" . The spatial filter guarantees 
that any SCB kernels that are detected are not associated with the flare ribbons. 
The resulting detections are shown in Figure [7j 



4. Physical Results of Algorithm Application 

The kernel tracking algorithm described in the previous section is applied to 
three different cases to train and refine the parameters. Table [l] outlines the 
basic characteristics of these datasets. From this initial application of the track- 
ing algorithm, several interesting results about flare and SCB kernels presented 
themselves. Section 4.1 discusses the derived physical properties of flare kernels 
and Section 4.2 characterizes the physical parameters of SCB kernels. Table [2] 
provides a summary of the results described in this section. 
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4.1. Properties of Flare Kernels 



We define a characteristic diffusion length L for flare kernels to maximize the 
length of trajectories formed while minimizing confusion between kernel tags (see 



Section 3.1.2). This parameter L corresponds to the maximum distance a kernel 
can travel between frames. A parameter of L = 4 pixels is found to be the best 
diffusion length for flare kernels. Thus, the maximum physical distance that a 
flare kernel can travel between frames translates to 3,200 km. Since L ~ |, where 
a is the characteristic distance between kernels, the diameter of the smallest 
resolvable kernel along the flare ribbon is approximately 6,400 km with the 
ISOON dataset. 

After flare kernels are associated between frames into trajectories, the tempo- 
ral integrity of flare kernels can be exposed. Flare kernels are temporally fairly 
robust; lasting on average ~ 70 minutes each. Generally, large flares take a few 
hours to evolve from the impulsive phase, through peak intensity, and return 
to pre- flare brightness. Also, the number of detectable kernels declines as the 
flare's intensity decays from its peak implying that there are fewer resolvable 
components in the flare ribbons as the flare evolves. This change in detectable 
kernels means that the majority of kernels cannot be tracked from pre-flare to 
post-flare suggesting a dynamic substructure to the flare ribbons when bright 
points appear and disappear as the flare erupts. 

Because flare kernels do not last for the entirety of the flare, individual kernels 
are not a good indication of overall flare behavior. However when taken in ag- 
gregate, flare kernels reproduce the overall topology of the flare. Figure [8] shows 
the intensities of all of the tracked flare kernels' plotted as points as a function 
of time. Integrating these intensities over each time step yields an aggregate 
intensity curve that reproduces the flare intensity curve shown in Figure [l]C. 

4.2. Properties of SCB Kernels 



Sequential chromospheric brightenings, although related to the erupting flare 
ribbons, are distinctly different than the flare kernels. Figure [9] plots the inten- 
sity of the entire population of SCBs versus time. It shows a relatively small 
time window when SCB occur before returning to pre-flare intensities. SCBs 
begin brightening about 30 minutes before flare peak and return to background 
intensity about 20 minutes after. In contrast, the flare intensity curve remains 
above pre-flare levels for several hours. 

Another result that is well demonstrated when plotting SCB kernel intensities 
in aggregate, Figure [9j is SCBs never have a peak intensity greater than their 
host flare. An SCB has a peak intensity 1.2 - 2.5 times above the background 
intensity level. In contrast, flares can brighten an order of magnitude or more 
above the pre-flare brightness. 

Individual SCBs are also much more fleeting than the flare ribbons. They 
increase to their peak brightness and return to background intensities in a much 
shorter timescale than the flare ribbons are observed to do. When isolated, these 
ephemeral phenomena have intensity curves that are impulsive with a sharp peak 
and then a return to background intensity in the span of about 12 minutes. 
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Table 2. A summary of results for flare and SCB kernels. 



Kernel 


Minimum Diameter 


Maximum Intensity 


Average Lifetime 


Motion 




(km) 


(above background) 


(minutes) 




Flare 


6,400 


10 


~ 70 


Directional 


SCB 


1,600 


1.5 


~ 12 


Random 



When the tracks of SCB kernels are examined, they do not show any progres- 
sive motion. The centroid of an SCB kernel randomly walks around within about 
6 pixels of its starting location for the duration of the trajectory. Although SCB's 
sequential nature of point brightening gives the appearance of a progressive 
traveling disturbance, the plasma beneath each brightening does not follow the 
disturbance and remains in the same location. This result confirms the findings of 
|Balasubramaniam et al. ( 2005 ) . The third row of images in Figure[7]are manually 
color-coded to reflect this traveling disturbance that is seen in an animation of 
the detections. The red group of detections matures quickly and dissipates while 
the yellow and blue groups take more time to develop. This suggests that there 
may be more than one population of SCBs. 



5. Discussion 

The method and results presented here suggest tracking flare kernels through the 
evolution of the erupting flare can characterize the evolving active region. We 
also demonstrate that a sum of the components of the flare ribbon reproduces 
the total intensity curve of the flare. Although it is not possible to say that 
any given kernel is tracking one specific flare loop, the flare kernels dissect the 
flare into its smallest visibly resolvable components in the ISOON Ha dataset. 
Since the number of detectable kernels declines as the flare's intensity decays, 
the overall intensity of the flare is related to the number of ribbon components 
discernible at any given time. 

The tracking algorithm also confirms the results of |Balasubramaniam et al. 



(2005) that sequential chromospheric brightenings are physically different from 
their associated flare ribbons. The tracked SCB kernels do not physically progress 
in one direction although the appearance of SCBs are sequential and propagate 
out from the flare ribbon. SCB kernels demonstrate an impulsive brightening 
without an exponential decay. Both of these physical characteristics of SCBs 
are distinctly different than the associated flare ribbons which do appear to 
physically move and have sustained intensity brightening lasting an hour or more. 
These differences imply a distinctly different physical mechanism is causing the 
sequential brightening as opposed to the flare ribbon brightening which confirms 



the assertion made by Pevtsov, Balasubramaniam, and Hock (2007). 



One major advantage of tracking flares and associated brightening with this 
method is that the locations and tracks of the flare and SCB kernels can be 
overlaid onto any dataset that is co-aligned with the original Ha images. A more 



SOLA: SCB_Techniques.tex; 8 August 2011; 1:29; p. 11 



Kirk et al. 



thorough and systematic examination of flares and SCBs is now possible allowing 
the user to develop statistics on any number of flares. Two ribbon flares without 
previously identified SCB's can also be analyzed with this automated routine to 
further investigate if SCB's are limited to moderate sized flares with associated 
CME's. To obtain a copy of this code, please contact the corresponding author 
with your request. 
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Figure 1. The two-ribbon flaring event from May 13, 2005. A shows an example of a calibrated 
Ha ISOON image with the region of interest (ROI) highlighted. The ISOON images are 2048 
x 2048 pixels with each pixel having a 1.1 arcsec resolution. B is the ROI after preprocessing. 
This is ROI covers 576 x 576 pixels corresponding to 2.12 X 10 11 km 2 on the solar surface. C 
shows intensity curves over the time period of interest: 13:00 - 20:02 UT. The ISOON Relative 
curve is an intensity curve defined by Equation ^ while the ISOON Peak curve is defined by 
the maximum value of the ROI in each time step; both curves are normalized. The GOES Soft 
(1.0 - 8.0 A) and Hard (0.5 - 4.0 A) curves are normalized x-ray intensity curves for reference 
from the GOES 14 satellite. 
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Figure 2. The two-ribbon flaring event from May 06, 2005. A shows an example of a calibrated 
Ha ISOON image with the region of interest (ROI) highlighted. The ISOON images are 2048 
x 2048 pixels with each pixel having a 1.1 arcsec resolution. B is the ROI after preprocessing. 
This is ROI covers 768 x 768 pixels corresponding to 3.77 X 10 11 km 2 on the solar surface. C 
shows intensity curves over the time period of interest: 12:28 - 22:47 UT. The ISOON Relative 
curve is an intensity curve defined by Equation ^ while the ISOON Peak curve is defined by 
the maximum value of the ROI in each time step; both curves are normalized. The GOES Soft 
(1.0 - 8.0 A) and Hard (0.5 - 4.0 A) curves are normalized x-ray intensity curves for reference 
from the GOES 14 satellite. 
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Figure 3. The two-ribbon flaring event from November 9,2004. A shows an example of a 
calibrated H« ISOON image with the region of interest (ROI) highlighted. The ISOON images 
are 2048 x 2048 pixels with each pixel having a 1.1 arcsec resolution. B is the ROI after 
preprocessing. This is ROI covers 600 x 600 pixels corresponding to 2.30 X 10 11 km 2 on the 
solar surface. C shows intensity curves over the time period of interest: 14:37 - 23:58 UT. The 
ISOON Relative curve is an intensity curve defined by Equation^while the ISOON Peak curve 
is defined by the maximum value of the ROI in each time step; both curves are normalized. 
The GOES Soft (1.0-8.0) and Hard (0.5-4.0 A) curves are normalized x-ray intensity curves 
for reference from the GOES 14 satellite. 
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Figure 4. A flowchart detailing the process in which flare ribbons and SCBs are identified and 
their physical parameters measured. The green parallelogram represents the database input. 
Blue trapezoids indicate a feature specific quality that must be determined and input. The 
gray diamonds are processes that have been defined by a collection of subroutines that are 
described in this paper and the blue shaded rectangles are individual routines. 
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Figure 5. A flowchart detailing the process in which flare and SCB kernels are identified. 
The green parallelogram represents the database input. Blue trapezoids indicate a feature 
specific quality that must be determined and input. The red diamond is a process that has a 
conditional result and the blue shaded rectangles are individual routines. 
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Figure 6. A time series of images showing the evolution of the flare ribbons for the May 13, 
2005 event. Each image shows the same flaring region ten minutes after the previous image. The 
flare peak occurs at 16:49 UT. The top row of images show the He? line center intensity. The 
second row of images show the flare after the ribbon isolation and enhancement as described 
in Section |3.2| The third row of images show the detected flare kernels after the image analysis 
algorithm is applied. The centroid of each flare kernel is marked by a black circle. Notice how 
the number of kernels change as the flare ribbon evolves. A comparison of the first and third 
row of images show the coverage of detections along the flare ribbons. 
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Figure 7. A time series of images showing the evolution of the SCBs for the May 13, 2005 
event. Each image shows the same flaring region three minutes after the previous image. The 
flare peak occurs at 16:49 UT. The top row of images show the He* line center intensity (same 
as the top row of Figure [61. The second row of images show a binary image of the same 
flaring region after applying the SCB isolation as described in Section |3.3| The third row of 
images show the detected SCB kernels after the image analysis algorithm is applied. Identified 
SCB's were visually separated into three groups to demonstrate how they appear to move 
when watching an animation of the flare. The red group propagates quickly while the blue and 
yellow groups develop more slowly. Notice how the number of kernels peak before the flare 
peak and rapidly declines as compared to the flare kernels. 
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Figure 8. A plot showing the May 13, 2005 measured intensities of each flare kernel plotted as 
points as a function of time. The solid blue line is the integrated kernel intensities at each time 
step. Overlaid for comparison in green is the normalized GOES 1.0 — 8.0 Aintensity curve. Here 
again, the vertical dashed line indicates the peak of the flare curve as shown in Figure[TJl!. The 
integrated flare kernels reproduces the shape of the flare curve determined by other methods. 
This indicates that the flare is well characterized by the flare kernels when taken in aggregate. 



SOLA: SCB_Techniques.tex; 8 August 2011; 1:29; p. 21 



Kirk et al. 



o o o o o 

o o o o o 

cc co ^ W O o 

^t 1 ^ ^ O 




Figure 9. A plot showing the May 13, 2005 measured intensities of each SCB curve plotted 
as points as a function of time. The solid yellow line is the integrated kernel intensities at each 
time step. Here again, the vertical dashed line indicates the peak of the flare curve as shown 
in Figure [l]C. Notice that the intensity of the integrated SCB curve returns to the continuum 
level much faster than the integrated flare kernel curve as shown in Figure [8] This indicates 
that the SCB kernels are measuring something distinctly different than the flare kernels. Also 
notice how the SCB intensities peak ~ 30 minutes before the flare intensity peaks. 
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